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We present a novel method to establish inner point contacts on hexagonal boron nitride (hBN) 
encapsulated graphene heterostructures with dimensions as small as 100 nm by pre-patterning the 
top-hBN in a separate step prior to dry-stacking. 2 and 4-terminal field effect measurements between 
different lead combinations are in qualitative agreement with an electrostatic model assuming point¬ 
like contacts. The measured contact resistances are 0.5-1.5 kfl per contact, which is quite low for 
such small contacts. By applying a perpendicular magnetic fields, an insulating behaviour in the 
quantum Hall regime was observed, as expected for inner contacts. The fabricated contacts are 
compatible with high mobility graphene structures and open up the field for the realization of 
several electron optical proposals. 


In recent years several experiments have shown, that 
ballistic graphene is an ideal platform for electron optical 
experiments. These experiments included the observa¬ 
tion of Fabry Perot resonances [1-3], snake states [4, 5], 
electron guiding [ 6 ], magnetic focusing [7, 8 ] or ballistic 
Josephson currents [9-12]. 

Up to now graphene encapsulated in hBN, which yields 
the highest quality for graphene on substrate, could only 
be accessed by top- or side-contacts from the edge of the 
device. [1, 13] However, a different type of contacts, in¬ 
ner point contacts (PCs) are required in order to realize 
several theoretical proposals on graphene such as e.g. the 
Veselago leasing [15] in single [16-18] layer graphene, val¬ 
ley [19, 20] or spin focussing [21] in graphene or for the in¬ 
vestigation of skipping and snake orbits of charge carriers 
at a pn-junction in combination with a strong magnetic 
field perpendicular to the graphene plane [22, 23]. 

In order to make PCs in the middle of the graphene 
sheet evaporated, sputtered or atomic layer deposited 
(ALD) dielectrics, such as MgO, Si 02 or AI 2 O 3 have 
been used so far [24]. However, these materials are in¬ 
ferior to the layered material hBN when it comes to the 
preservation of the graphene quality [25]. It is possible 
to establish PCs on graphene using a STM tip where also 
the position of the contact is changeable. However, do¬ 
ing this at low temperatures, involving several PCs at 
the same time is extremely challenging. 

Here we present a novel method to establish PCs to 
graphene encapsulated in hBN by pre-patterning the h- 
BN flake before the stacking process. This allows to 
access the graphene at arbitrary position with contacts 
smaller than 100 nm in diameter. The method is compat¬ 
ible with clean graphene fabrication, since the graphene 
transport channel will not be in contact with any resists 
or solvents during fabrication [ 1 , 26-28]. We extract the 
graphene quality by comparing measured 2 and 4 termi¬ 


nal resistance values with a simple model. Furthermore 
we show localization of the edge-states around the PCs 
in a magnetic field, expected for a proper inner contact 
in the quantum Hall regime. 



FIG. 1. Fabrication of hBN-graphene-hBN hetero¬ 
structures with top-PCs. a, False color SEM image of 
the top-hBN on Si02 substrate after drilling the holes with 
the Ga-FIB. Inset: Glose-up of a single hole having a diameter 
of ~100 nm. b, Optical image of a final stack. The different 
layers are indicated in red (bottom hBN), black (graphene) 
and blue (top hBN). The holes are indicated with arrows, c, 
False color SEM image of the final stack (blue) with Pd con¬ 
tacts (yellow) overlapping the drilled holes, d, Schematic of 
the cross-section as indicated in (c) with pink arrows and the 
dashed line. 

To produce PCs compatible with high mobility 
graphene we have fabricated four equally spaced holes 
in the top hBN layer of hBN-graphene-hBN hetero¬ 
structures. In order not to damage the graphene, the 
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top-hBN layer (^10-20 nm thick) was pre-patterned on 
a separate Si 02 wafer using a gallium based focused ion 
beam (Ga-FIB, see details in supporting) prior to the 
dry-stacking of the heterostructure [1, 29]. In contrast to 
establishing the holes with conventional e-beam lithogra¬ 
phy and subsequent etching, the drilled holes are better 
defined in shape and the diameter can be adjusted more 
reliably. In our samples we use a hole-diameter of ap¬ 
proximately 100 nm as it is shown in Fig. la. Before 
picking-up the top-hBN from the Si02 wafer, it is briefly 
exposed to a CHF 3 /O 2 plasma treatment as without the 
plasma the flakes are pinned to the Si 02 and pick-up is 
not possible. 

The top-hBN flake is transferred from the Si02 sub¬ 
strate to a 1 jam thick poly-propylene-carbonate (PPG) 
polymer by spinning the PPG directly onto the Si02 chip 
containing the top-hBN with the holes and then peeling 
the polymer off. The remaining assembly procedure of 
the hBN-graphene-hBN heterostructure follows the dry¬ 
stacking approach proposed by L. Wang et al [1]. Since 
only the top side of hBN comes in contact with the poly¬ 
mer, the method preserves the clean fabrication of dry¬ 
stacking graphene. The bottom hBN flake is exfoliated 
to the target wafer directly, with typical thickness of 20- 
30 nm. 

We use a strongly doped Si++/Si02 wafer with a 300 
nm thick oxide to gate our devices. The final stack was 
annealed in forming gas at T=300 °G for 3 hours in or¬ 
der to reduce strain and maximize the areas without 
bubbles [30]. A contrast adjusted optical image of the 
annealed stack is shown in Fig. lb. The 100 nm thick 
palladium (Pd) contacts are established using standard 
e-beam lithography and e-gun evaporation. A false color 
SEM image of the contact area is shown in Fig. Ic. A 
schematic of the cross-section of the stack with contacts, 
as indicated with the dashed pink line in Fig. Ic, is shown 
in Fig. Id. Further details of the fabrication are given in 
the Supporting Material. 

In total 4 different samples were produced all showing 
a similar behaviour. The measurements were performed 
at cryogenic temperatures (1.5-4 K) using standard low- 
frequency lock-in technique. 

In the following, first the contact resistance and the 
field-effect measurements at zero magnetic field are dis¬ 
cussed. In the second part, the behaviour of the devices 
at high magnetic field perpendicular to the graphene 
plane is presented. For the calculation of the charge 
carrier mobility, a model fitting the device geometry is 
introduced. The nomenclature of the various differential 
resistances is given as Rij^nm = dVnm/dlij^ with lij flows 
from PG i ^ j and the voltage Vnm measured as the dif¬ 
ference between PG n and m, with i, j, n, m G 1 — 4. 

Figure 2a shows Rij^nm foi* all possible 4-terminal con¬ 
figurations. Out of the six possible configurations, only 
three are independent: measurements where current- and 
voltage-probes are inverted are identical as expected from 


the Onsager relations [31]. The resistance traces show a 
sharp maximum around zero doping, corresponding to 
the charge neutrality point (GNP) of graphene. 

For rectangular graphene devices, where the current 
density within the graphene sheet is constant, the mo¬ 
bility ja can be deduced by measuring the field effect of 
the longitudinal resistance Rxx{yBG)i taking the length 
and width of the Hall bar into account. For PGs, which 
are situated in the middle of the graphene sheet, a dif¬ 
ferent formula has to be used since the current density 
within the graphene sheet varies. Here we introduce a 
model to extract the sheet conductivity (a) from the 4- 
terminal measurement of the resistance, assuming an infi¬ 
nite graphene sheet with a constant a. The four contacts 
are at positions Tx with x = i, j, n, m and diffusive trans¬ 
port. Starting from a single PG at position ri, the current 
spreads isotropically into the graphene. This leads to a 
current density of j(r) = //(27r|r —rij) at distance |r — rij 
away from the PG, where / is the current. According to 
j(r) = E(r)cr, the electric field, E(r) ^ l/|r — rij, leading 
to an electrostatic potential H(r) ^ In (|r — rij). Assum¬ 
ing a current flow between two PGs from i ^ j, the 
potential at position r is obtained by the superposition 
principle: 

nr) = -^*ln(|r-ri|) + ^*ln(|r-rj|) + C, (1) 

where G is an integration constant. In the 4-terminal 
measurement only the voltage difference between the two 
leads at position n and m {Vnm = V (i*n) —E(rjn)) is mea¬ 
sured. For simplicity we assume that the voltage probes 
do not influence the electric field pattern in graphene as 
shown in Fig. 3a and b. This leads to 


Vnm = * 111 

27rcr 


/ |rn-rj| |rm-ri| \ 

V|rn-ri| Ir^-rjiy ’ 


( 2 ) 


where the conductivity a can now be extracted from the 
measurement. 

Using Rij^nm O' can be deduced. The mobility, ja of the 
graphene was then extracted using the linear dependence 
of a{n) on the carrier density, n {a = neji). The density 
was calculated from the back gate voltage using a 
parallel plate capacitor model. The hole and electron 
doped region revealed mobilities of jih r\j 35’000cmV(Vs) 
and jie ^25’000 cm^/(Vs) respectively. This is in 
good agreement with the conversion extracted from 
the evolution of the filling factors in a quantum Hall 
experiment (QHE) (see supplementary information). 
An alternative way to determine the mobility is by the 
onset of the Shubnikov de Haas (SdH) oscillations in the 
QHE measurement. According to jiB = uJcT = 1, where 
ujc is the cyclotron frequency and r is the scattering 
time, the charge carriers can complete a full cyclotron 
orbit before being scattered. With the SdH oscillations 
starting at ^0.4-0.5T, a corresponding mobility of 



3 



Vbg (V) V,, (V) 


FIG. 2. 4-terminal measurement of the graphene re¬ 
sistance between different combinations of PCs. a, 

The original 4-terminal resistance Rij^nm measured between 
the six different lead combinations. Measurements with in¬ 
verse voltage and current probes are identical, b, Resistance 
of the same data presented in (a) multiplied by a factor de¬ 
pending on the measurement geometry as given by Eq. 2. In 
the non-local measurement Ri 2 , 34 , the voltage V — Vn — Vm 
was negative in the n-doped region. 

11 ^20’000-25’000cm^/(Vs) is extracted which is in 
good agreement with the values deduced from the field 
effect measurements. Confirmation that a single layer 
graphene (SLG) is encapsulated in hBN was given by 
the observed sequence of filling factors in magnetic 
field (see supplementary information) and by Raman 
spectroscopy (see supplementary information). 

From Eq. 2 it follows that in case of a homogeneous a, 
all 4-terminal resistance measurements can be renormal¬ 
ized according to 


much closer to each other than the unsealed ones, which 
confirms that our theoretical model is realistic. 

The deviations from the ideal case can be related to 
the boundary conditions assumed for the model. The 
most significant deviations from the ideal model are 
probably i) the finite dimensions of the metallic PCs 
and ii) the finite size of the graphene sheet, which both 
change the electric field pattern. Besides that, the sheet 
conductivity does not seem to be fully uniform within 
the sample as can be seen by the slight shift of the 
charge neutrality point between several measurements. 
The charge neutrality points are at I/=-2.6V, -2.8 V 
and -4.8 V, respectively. Moreover, a non-uniformity 
of the doping profile arises also from the screening of 
the top contacts. This results in different lever arms 
of the back-gate for regions covered and not-covered 
by the electrodes. Finally, Eq. 3 was derived assuming 
a completely diffusive sample. However, the charge 
carriers in the sample are most likely in an intermediate 
regime between the diffusive and the ballistic regime. 
Using the Drude formula, the scattering mean free 
path can reach 1 jim at Vbg = 30 V, which is in the 
same order as the contact distance a = 2.2 jim. In this 
intermediate regime it can occur that the voltage drop 
over a larger, but clean (ballistic) area is lower compared 
to a smaller, but dirty (diffusive) area. Moreover, 
for ballistic trajectories the probability of arriving at 
a contact, which is farther away can be higher. This 
picture explains the negative resistance at certain doping 
values observed in the non-local measurement indicated 
with an arrow in Fig. 2c. For all the configurations 
where the voltage-probes are (at least partially) within 
the current path, the bias voltage will be dominant and 
consequently no negative signal can be seen. 


Rij,nm — Rij,nm 


/ |rn-rj| |rm-ri| 
V|rn-ri| Ir^-rjl 


-1 


(3) 


With all four contacts at equidistant spacing the 
logarithm in Eq. 3 simplifies to ln(4), ln(3) or ln(3/4) 
depending on the measurement configuration. R the 
model with all the assumptions is valid, all Rij^nm 
should be equal and given by l/27rcr. The difference 
between the original and renormalized values can be 
seen in Fig. 2a and 2b, respectively. 


The renormalized values Rij^nm in Fig- 2b are not ex¬ 
actly overlapping, as would be expected for a perfect sys¬ 
tem given in Eq. 3. However, one can see that the non¬ 
local measurements (voltage probes outside the current 
path), shown in green, which were in the original data 
smaller by a factor of 7.5 (8) from the blue (red) curve, 
deviates now only by a factor of 2 (1.6) after renormal¬ 
ization. On the other hand, the local measurements, the 
blue and red curves are in rather good agreement be¬ 
fore and after renormalization. Overall, the rescaled are 


In order to extract the contact resistance Rc which 
arises between the metal leads and the graphene, Rc 
(interface resistance), we turn to two terminal measure¬ 
ments. To calculate the contact resistance we measure 
the 2-terminal resistance R 2 T between the outer contacts 
(1,4) as sketched in Fig. 3a. Then the contact resistance 
can be calculated according to: 

R 2 T = ‘^Rc + ^ • Rat^ (4) 

where A • R^t = A • i?i 4,23 is the intrinsic graphene 
resistance including a geometry factor A (which will 
be evaluated in the following) and R 2 T = Ria,iA' Here 
we assumed the same contact resistance for the two 
contacts. Due to the higher electric field near the source 
and drain contacts, shown in Fig. 3a, the voltage changes 
faster near them. This can be seen in Fig. 3b, where 
the potential along the line connecting the contacts is 
plotted. The resistance coming from the non-linearity 
of the potential near the contacts is called spreading 
resistance and leads to a potential difference that is 
marked as Vsp in Fig. 3b [32]. It is of the same origin 
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FIG. 3. Contact resistance extracted from 2- and 4- 
terminal measurements, a, The 2-terminal measurement 
setup and a sketch of the corresponding electric field pro¬ 
file. The 2-terminal measurement is performed between the 
two outer-most PCs which have a distance of 3a from each 
other, whereas the reference four terminal measurement be¬ 
tween contacts 2 and 3. b, The red, dashed line represents 
the linear interpolation of V23 compared to the actual electro¬ 
static potential (black, solid line) along the line connecting 
the contacts. The resulting potential difference at distance 
d/2 from source- and drain-center is indicated with Vsp. c, 2- 
terminal resistance between 1 and 4 is shown in red, whereas 
the 4-terminal resistance i^i 4,23 is shown in blue. Extracted 
contact resistance of the configuration given in (a) using Eq. 4 
and 6 


as Maxwell’s resistance that occurs in metallic point 
contacts [33]. 


To calculate the geometrical factor A, we use Eq. 2, and 
apply it together with Eq. 4 to the situation sketched in 
Fig. 3a. Because the potential is singular at position ri 
and Tj we introduce a cut-off for the potential at (i/2 away 
from the singularity, where d is the diameter of the con¬ 
tact. Physically speaking this accounts for the equipo- 
tential within the metallic PCs with a finite dimension. 
Doing so, both terms in the numerator and denominator 
of Eq. 2 become (3a — (i/2) and d/2 respectively, where 
a is the distance in between two neighbouring contacts. 
This leads to 


Using R^t = iii 4,23 where iii 4,23 = /n(4)/(27rcr) the ge¬ 
ometrical factor becomes: 


A = In 




1 

In (4) ’ 


(6) 


In Fig. 3c the extracted contact resistance is shown using 
a geometry factor of ^4=7.04 (a=2.2 /im and d=100 nm). 
The contact resistance of different contact configurations 


and different devices is of the order of Rc=0.5-1.5 kU 
at high doping. This value is quite remarkable for PCs 
of only 100 nm in diameter since as well top-contacts 
with significantly larger areas (in the order of jam^) have 
resistances in the kU range. Moreover, the presented 
model shows that the contact resistance at high doping is 
given roughly by the 2-terminal resistance at high doping 
(UBG==t30 V). In this case graphene becomes very con¬ 
ductive and the voltage-drop over the graphene is mini¬ 
mal, therefore R 2 T ^ 




B(T) 


FIG. 4. 2-terminal conductance as a function of back- 
gate and magnetic field, a The conductance between 
neighbouring PGs (distance of 1 /xm) becomes zero at high 
magnetic fields. The black area shows the threshold for 
G < 0.025 e^/h (R >1 MQ). A line-trace at R = 15 T is 
given on the right-hand side of the color-plot, b, Numerical 
derivative calculated with respect to Vbg of the measurement 
given in (a) reveals the constant-conductance regions more 
pronouncedly. 

To further characterize the PCs, we have applied 
a magnetic field perpendicular to the graphene sheet, 
which forces the charge carriers to move along cyclotron 
orbits. With a sufficiently high magnetic field the de¬ 
vice is driven into the quantum Hall regime, into a state, 
where the bulk of the sample is insulating, since charge 
carriers will be localized either around the PCs or along 
the edges of the sample, which decouples the PCs from 
each other and the edge of the sample. In the case of 
a homogenously doped and gated device one would ex¬ 
pect complete insulation as soon as the cyclotron orbit 
and the magnetic length are smaller than the distance 
between the PCs. The latter is the case for magnetic 
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fields in the few hundred mT range. Moreover > 1 
is required. 

The conductance as a function of back-gate and mag¬ 
netic field of a sample with a hole to hole spacing of 
a=l /rm and a mobility of /i =15’000 cm^/(Vs) is shown 
in Fig. 4. The magnetic field dependence of the device 
discussed before can be found in the Supporting Infor¬ 
mation. The black region in Fig. 4a shows values of 
Go < 0.025 e^/h {R >1 Mf7). It can be seen, that by 
the application of magnetic field the sample becomes in¬ 
sulating, however, the fields required are much higher 
than expected. From simple considerations above this 
should happen around 0.5 — IT. 

Inhomogenous doping distribution in the sample can 
result in coupling of the contacts. For complete insula¬ 
tion of the device, the graphene has to be simultaneously 
insulating in the whole region between the contacts, since 
on the border of regions with different filling factors edge 
currents will flow. Tunneling to these edge states can give 
short-cut currents between the contacts. As already men¬ 
tioned in the text, doping inhomogenities exist within our 
sample, as can be seen from the shift of the charge neu¬ 
trality point between different 2T measurements. Fur¬ 
thermore, locally the back-gate can be screened by the 
top contacts. However, our estimates have shown, that 
the screening changes the gate efficiency by less than 4% 
far away from the CNF, where the quantum capacitance 
is high. Only close to the CNF, where the quantum ca¬ 
pacitance is small, does screening of the contacts increase 
to 20%. Furthermore, an additional offset potential may 
emerge in the regions of the top contacts due the forma¬ 
tion of a contact potential between the palladium con¬ 
tacts and h-BN. We emphasize, that substantial part of 
the voltage drops in the region close to the contact. The 
combination of all these effects can cause local differences 
in the filling factor, which can account for the observed 
high threshold fields. 

In future devices insulation of at lower field can 
be achieved by choosing devices with smaller doping 
inhomogenities, which can come from bubbles present in 
the stacks. Moreover the inhomogenous screening of the 
top contacts or offset potentials can be circumvented by 
careful design, in which the fiake would be fully covered 
with a metallic plane to achieve a homogeneous doping 
situation. 

We have shown a new method to establish inner point 
contacts with dimensions of 100 nm in a hBN-graphene- 
hBN heterostructure. A simple model has been intro¬ 
duced which qualitatively explains our 2 and 4-terminal 
gate dependent conductance measurements. Surprisingly 
low contact resistance Rc=0.5-1.5 kQ have been found 
despite the small FC size. Magnetic field measurements 
showed that the inner contacts are decoupled from the 
edge and from each other at high magnetic fields. 

The presented technique is compatible with high- 


quality encapsulated graphene, since the hBN fiake 
is patterned prior to the stacking and therefore the 
graphene remains clean. With further optimization 
one can expect devices with mobilities around 100 000 
cm^/(Vs). The technique also holds the potential to fur¬ 
ther decrease the contact size, since with the Ga-FIB hole 
diameters below d=20 nm are possible. 

The point contacts introduced here give the possibility 
to complement side and top contacts in complex devices 
and could be the potential milestone towards realizing 
novel concepts like lensing or measurements of caustics 
in p-n junctions. 
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Fabrication 
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FIG. SI. Fabrication of the point contacts on a hBN-graphene-hBN heterostructure (not to scale), a, Drilling 
the holes into the top-hBN using a Ga-FIB. b, Due to the drilling, Si02, hBN and Ga is deposited on the side of the hole and 
in the near proximity of the hole on top of the hBN (indicated with purple dots). Exposing the top-hBN flake to a GHF 3 /O 2 
plasma avoids pinning of flake to the Si02 substrate, c, The top-hBN flake with the holes is removed from the Si02 support 
by spin-coating PPG on top of the wafer and then peeling it gently of. The remaining assembly of the hBN-graphene-hBN 
stack is shown in the figures (d)- (e) where the black arrow indicate the pick-up of graphene and the release of the half-stack 
on a bottom hBN respectively. After complete assembly of the stack, the PPG is dissolved in Ghloroform. f, The palladium 
contacts are assembled using standard e-beam lithography and e-gun evaporation. 


In order to drill holes into the top-hBN with a gallium based focused ion beam (Ga-FIB) we use an acceleration 
voltage of 30 keV and the smallest possible current (1.1 pA) in order to obtain highest resolution. The hBN to be 
patterned was exfoliated on a Si++/Si 02 substrate with a 315 nm thick oxide, using the scotch-tape technique. The 
chips were previously carefully cleaned using Piranha solution (98% H 2 SO 4 and 30% H 2 O 2 in a ratio of 3:1) since it 
is the bottom face of the hBN which will later on contact the graphene. Once an ideal hBN flake (thickness ^10-30 
nm) is identified by optical microscopy, the Ga-FIB is used to drill several holes into the flake (with diameter d ^100 
nm and a equidistant spacing of 1 - 2.2 fim) as sketched in fig Sla. 

Before picking-up the hBN from the Si 02 wafer, it is briefly exposed to a CHF3/O2 plasma (40 sccm/4sccm, 60 
mTorr, 60 W, 15 s) as shown in fig. Sib. It turned out that without exposing the hBN flakes to the plasma, the hBN 
flakes could not be picked-up from the Si 02 substrate. A possible explanation might be that during the drilling of 
the holes with the Ga-FIB, Si 02 from the wafer is sputtered on the side of the holes which pins the flake to the wafer. 
The CHF3/O2 plasma removes this layer and allows therefore a successful pick-up of the flake from the Si 02 chip. 
To pick-up the top-hBN, the Si 02 chip is spin-coated with ^1 /im of poly-propylene-carbonate (PPG) and baked 
at 80 °C for 5 minutes. By peeling-off the PPG gently from the substrate (fig. Sic), all hBN flakes are transferred 
from the Si 02 onto the PPG polymer. Peeling-off the PPG without breaking the hBN flakes works best when slowly 
releasing the PPG at a low angle from the Si 02 chip (drilled flakes are more likely to break). The PPG with the hBN 
flake is then placed on a home-made stamp of ^0.5 mm PDMS. 

The remaining assembly procedure of the hBN-graphene-hBN heterostructure follows the dry-stacking approach pro¬ 
posed by L. Wang et al as shown in fig. Sld-e. [1]. The 100 nm thick palladium (Pd) contacts are established using 
standard e-beam lithography and e-gun evaporation. A cross-sectional sketch of the final device is shown in fig. Slf. 





































Low field measurements 



FIG. S2. 4-terminal measurements as a function of back-gate and magnetic field, a, Derivative djdVBG ^ 12,34 of 
a non-local measurement. The sequence of the filling-factors (red, dashed lines), z/ = di 2, ± 6, di 10... is in agreement with 
the ones expected for SLG. b, 4-terminal measurement of the resistance i?i 4,23 which is more comparable to a classical QHE. 
A line-cut at T is given at the right-hand side of the color-plot. The corresponding filing-factors are indicated. 

In a Hall bar configuration, a clear distinction between longitudinal- {Rxx) and Hall-resistance {Rxy) can be made. 
As long as the Fermi energy is in between two Landau levels (LL) Rxx = 0 while Rxy = {hjne^). However, in our 
sample the situation is more complex due to the absence of a graphene edge which directly couples to the contacts. 
As all four contacts are situated in a row, a separation between longitudinal- and Hall-resistance is impossible. 
Therefore, the 4-terminal resistance i?i 4,23 5 shown in fig. S2b, is more complicated to interpret and goes beyond the 
scope of this studies. The filling factors have been assigned based on a capacitance model. 

The evolution of the filling factors with varying back-gate and magnetic field was determined using a non-local 
measurement as it revealed more pronounced features in the fan-plot. The back-gate voltage to density conversion 
(n = oVbGi where a is the lever-arm of the back-gate) was extracted from the evolution of the filling factors 
{Yeg = Bevl{ah)^ where v is the filling factor of the LL and h is the Planck constant) in the 4-terminal, non-local 
measurements shown in fig. S2a. The evolution of the filling factors (z^ = ± 2, ± 6, ± 10...) shown in fig. S2b are in 
good agreement with the expected sequence for single layer graphene, indicated with the red, dashed lines. 
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FIG. S3. Raman spectra of graphene encapsulated in hBN. The exact number of graphene layers encapsulated in 
hBN was determined using Raman spectroscopy. The FWHM of the 2D-peak at roughly 2690 cm“^ was fitted with a single 
Lorenzian (blue curve in the inset). The extracted FWHM of 18.3 cm“^ fits best to the value expected for single layer graphene. 


Fitting the full width at half maximum (FWHM) of the 2D-peak in the Raman spectra is an alternative way to 
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determine the exact number of graphene layers present in the device, shown in fig. S3. [2] We extracted a FWHM 
of 18.3 cm“^ using a single Lorentzian to fit our data which fits much better to the value of single layer graphene 
(FWHM=27.5±3.8 cm“^) as compared to bilayer graphene (FWHM=51.7.5±1.7 cm“^). 

Having more knowledge about the evolution of the filling factors with back-gate and magnetic field, the various 
filling factors could be assigned to the peaks and valleys in fig. S2b. A line-trace at B=b T is given on the right-hand 
side of fig. S2. 


Thermal cycling 



FIG. S4. Thermal cycling of point-contact devices The devices did not change much during thermal cycling between the 
first and second cool-down as can be seen on the two measurements of the 4-terminal differential resistance. Thermal cycling 
was required to transfer the devices from the first measurement setup (B <8 T) into a setup with high magnetic fields (B <15 
T). 

The identical 4-terminal measurement was performed after the first- and second cool-down to obtain some knowledge 
about a possible degradation due to thermal cycling. As can be seen in fig. S4, almost no change could be observed 
between the two measurements. 
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High field measurements 



FIG. S5. 2-terminal conductance as a function of back-gate and magnetic field, a The conductance between 
neighbouring PCs (distance of 2.2 /xm) becomes insulating at high magnetic fields. A line-trace at B=15 T is given on the 
right-hand side of the color-plot, b, Numerical derivative (d/dVeG Go(e^//i)) of the measurement given in (a). 


High-field measurements of another sample with a distance of 2.2 jam in between the PCs is shown in fig. S5. It 
can be seen that the sample becomes insulating between neighbouring contacts at comparable fields as the sample 
shown in the main article. 

The black region in fig. S5a shows the threshold for Go < 0.025 e^/h {R >1 MQ). 


Christian.Schoenenberger@unibas.ch 

[1] L. Wang, I. Meric, P. Y. Huang, Q. Gao, Y. Gao, H. Tran, T. Taniguchi, K. Watanabe, L. M. Campos, D. A. Muller, 
J. Guo, P. Kim, J. Hone, K. L. Shepard, and C. R. Dean. One-dimensional electrical contact to a two-dimensional material. 
Science, 342(6158) :614, November 2013. 

[2] Yufeng Hao, Yingying Wang, Lei Wang, Zhenhua Ni, Ziqian Wang, Rui Wang, Ghee Keong Koo, Zexiang Shen, and John 
T. L. Thong. Probing layer number and stacking order of few-layer graphene by raman spectroscopy. Small, 6:195, 2010. 



















